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Abstract 

Following the general form for LQG optimal compensators for flexible 
structures with collocated rate sensors we develop an explicit compensator 
transfer function for the NASA LaRC CSI Evolutionary model in the form: 

y(i'co) = gi(iyB*(-M b co 2 + 7’(i(o) + iy(oB u B*)~ l B u 

where T(jco) is a 48x48 positive definite matrix whose derivation is the 
main result of this report. The undamped mode frequencies can be expressed 
in terms of T(i(o) as the zeros of 

Det \-Oi 2 M b + T(/(0)| 

while “clamped-clamped” modes of the structure (with all nodes clamped) 
are the poles. 


1. Introduction 

In this paper we present an explicit compensator transfer function for the NASA 
LaRC Evolutionary Model [1], using the Continuum Model developed in [2], In particular 
the notation follows closely that in [2]. The compensator is obtained upon specialization 
of the general development in [3], 
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2. Compensator Transfer Function 

The optimal compensator transfer function is given by (see [3]): 

-£= 'id a ld r B*(p 2 M + A + pyBB*)~ ] B . 

The main step is to calculate 

V| f{p) = B*(p 2 M + A + pyBB*)~ l B . 


We shall consider only the Continuum Model in [2] (case 3) in which the main bus, the 
tower and the appendages are flexible but the antenna is lumped. 

Let 


y(p)v = u 

B*(p 2 M + A + pyBB*)~ l Bv = u . 


Let 


( p 2 M + A + pyBB*)~ X B\ 


f 

b 


Then 


p 2 M 0 f + A 0 f = 0 


where / is also subject to the “linkage conditions” (see 12]), 

p 2 M b b + A b f + pyB u B*b = B u \ 
and 


( 2 . 1 ) 


( 2 . 2 ) 


u = B*b . 


We shall now specialize to p = i(0, to real. 
To solve (2.1), we let 


Zl 


z = 


28 
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where 


Let 


if p = /O), 


z i = /( 0, 0, 0) 
z 2 = /,( 0, 0, 0) 
z 3 = fy(s 2, 0-, 0) 
z 4 = / y (j2. 0+, 0) 
z 5 = f z (s T , 0, 0+) 
z 6 = /*(J4+. 0, 0) 
z 7 = f y {s 5 , 0— , 0) 
z 8 = fy(s s, 0+, 0) 


0 




^2 (^0 + P 2 ^ o ) ^ 2 '^! 


A~p 2 ) = jtf(co) 


in the notation of [2], From now on, let p = i'oj. Let 


J(w)s 

e 


Pn(s) P l2 (s) 
Puis) P 22 {s) 


Ps 2,1 l( 5 ) Ps j,12 i s ) 


and as before in [2], Then 


/(*, 0, 0) 

A (co )s 


f x (s, 0, 0) 

= e 

*2 


0 < 5 < 52 
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and 


/(i 2, 0, 0) 


*1 

/.(i2-, 0, 0) 

= e 

*2 


f{s 2 , S, 0 ) 
fy(S 2 , S, 0) 


e A S t ^ s 




f(s 2 ,0,0) 

*3 


f(s 2 , 0,0) 

Z4 


< 5 < 0 


o < j < e, . 


Let us now display the values at the nodes only. Then going along the main bus: by 
Linkage Condition (2) (in [2]): 


f x (s 2 +, 0, 0) - f x (s 2 ~, 0, 0) + A 2 1 A 2 ' Si (zj -z 4 ) . 


Hence 

f( s T< 0 » 0 ) 

f x ( s T~' 0, 0) 
and 

/(j 4 , 0, 0) 

f x {s 4 ~, 0, 0) 

By Linkage Condition (1) (in [2]): 


e 1 


/to, o.O) 

f x (s 2 +, 0, 0) 


^(co)(j 4 -J r ) 


/(5 r ,0,0) 

f x ( s T+> 


f x (s T +, 0, 0) - f x {s T -, 0,0) - A~ 2 l {L l T f{s T ,0,0) - A 2 J z s ) 


/(is, 0. 0) 

/.(is-, 0.0) 


4(a>)(j,-* 4 ) 

e 


/(i 4 ,0,0) 

26 


and finally 

/(is, 0, 0) 
/.(is+, 0, 0) 


/(L, 0, 0) 

i(co)(L-j 5 ) 

/,(L, 0, 0) 

= e 
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where by Linkage Condition (3) (in [2]): 


f x (s$+, 0, 0) = f x (s$—, 0, 0) + '42,s j ( z 7 — *8) • 

Next: 


f( s T< ^t) 

(CO)t, r 

— /p 1 

/(*t,0,0) 

f z ( s T’ Lj) 

— t 

z 5 

f( s 5. ~^2. 0) 


z) /(* 5 . 0 , 0 ) 

~^2i 0) 

— 6 

z 7 

y(- y 5 » +^2t o) 

/? J 

/(5 S ,0,0) 

/ z (j 5 , +e 2 , 0) 

— C 

Zg 


Hence we can calculate /(■,•,•) as 

/ = £(co)z . 

In particular we can calculate b in terms of z. Let 

b = L fc (co)z 

where L b { to) is a (matrix) function of co. Let b = col. (b\, .... b%) and 

L b = {L iy } , each L iy being 6x6 . 

Then 

b\ = zi ; L\\ = / , L l( = 0 , i * 1 

bi = /(*2, -*i, 0) 

= ( _ ^l)/( J 2> 0. 0) + ^l) z 3 

= 1 ( — ll( J 2) z l + 2 C — ^1 ) z 3 • 



Hence 

^2\ = ^5 2 ,l 1 ( — 1 1 (^ 2 ) 

^22 = ^ 2 ,11 (-*1^120*2) 

^23 = Ps lf 12(“^l) 

^3 “ f( s 2* + ^1* 0) = ^5 2 ,n ( + ^l)(^ > ll( J 2) z l + P\l{ s l) 2 l) + P S 2,12 (^1 ) z 4 

^31 = /^.iitfiFiita) 

^32 * P Sli l 12( J 2) 

^33 = 0 

^34 = P 5 2> 1 2(^l) 

t>4 = /(^r* 0, ^t) 

= ^ 5 r ,i i(^r)[^n( J r - * 2 ) 0 Pnte)*i +^120*2)22) 

+ P 12(^7 “ S 2)(P 2l(^2) z l + P 2 2 (^2)^2 + ^2 ^2,5 2 ( z 3 ” z 4))] 
+ ^ 5 r> 12 (^r ) z 5 

^41 = ^S T$ \ 1 1 l( J 7* - s j)P\\( s 2) ^ P\2( S T ~ s 2)P\2( s 2)] 

L 42 - Ps Tt \l(lv)[P\d s T - S 2)Pl2(s 2 ) + P\2(S T ~S2)P2l(S2)] 

^43 = Ps Tt \\ih')(P 12(^7 “ ^2)^2 ^2,S 2 ) 

L 44 = ~Ps ft \ \(Lt)P\ 2 ( s T ~ s 2)^2 A 2 f s 2 

^45 = ^5 r ,12(^r) 
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b s = f(s 4t 0, 0) 


^51 “ ^11(^4 “ ^ 7 *) ^11 ( S T “ * y 2)^ > n(* y 2) + ^ 1 \( s 4 “ s T)Pl2( s T ” 5 2)^2l(^2) 

+ P 1 l( s 4 “ s t)P 2\( s T ~ s 2 )P 1 1(^2) + P n( s A - s t)P 22 ( S T ~ s 2)Pl\( s 2) 

— ^12(^4 “ s t)P 22( s T ~ s 2)^2 j{Pn( s T~ s 2)P ll( 5 2 )+P 12(- y r“* y 2)/ 5 2l(‘ y 2)} 

^52 = P 11(^4 ~ s t)P u( s t “ J 2)^ 12(^2) + ^ nfar s 2 )P 22(^2) 

+ P 12( J 4 ” s t)P 21 ( S T “ 5 2)^ 12(^2) + P 12(^4 “ $t)P 22( s T ~ S 2 )P 22(^2) 

- P 12(^4 — s t ) P 2 i( s t ~~ **2)^2 ^1 ,r{^ ii(* y r"'* y 2)/ > 12(^2) + P \ i ( s t ”**2)^2 2 (^2)} 

^53 = ^ll(^4 - ^r)^ > 12(* y r ” ^2)^2 ^2,S 2 + ^12(^4 - *y7')^ > 22(- y r ” ^2)^2 ^2 ,S 2 
- ^ 12(^4 “ s T ) P 22 ( J r _ **2)^2 ^1 ,r{^ ” s 2)^2 A 2 t s 2 } 

^54 = -P\l(U ~ S T )P n (s T - s l) A 2 X Al t S 7 

+ ^ 12(^4 “ P22( S T “ J 2)^2 ^1 ,r {^ 12(^7 ~ ^ 2)^2 ^2 ,s 2 } 

£55 = ^ 12(^4 - 2li s T ~ s 2)^2 A 2 j 

be = f ( s S > 2» 0) = ^5 5f i 1 (~^2)( P \ i ( s 5 “ ^4)^5 + ^12(^5 - J 4) z fi) + Ps 5 , 1 1 ( ^2 ) z 7 


Ui - 

^ 5 5 ,1 1 (”^2)^ 11(^4 ’ 

•^5)^5 1 » 

1 < 5 

^66 = 

P s 5t 11 ir ^ t)P 12(^4 

- *5) 


^67 = 

Ps 5 m ^ 2 ) 



^7 

= f ( ss * + ^ 2 > 0 ) 



= 

P s 5,1 1 (^2)^ 11(^5 “ 

‘^4 )^5 i > 

vn 

VI 

^76 = 

P S 5,1 1 (^2 )P 12(^5 - 

- *4) 


L77 = 

0 



^78 = 

^55,12(^2) 
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h = f{L, 0, 0) 


^81 = 1 1 — • y 5)^ > 1 xC-^S — • s 4) + ^ 12(^- — -^5)^* 2l(- y 5 ~ -^4)1^5 / » 1 ^ 5 

1*6 = P 1 1 (^ - s s)P n( 5 5 - ^ 4 ) + P n(L - ^ 5 )P 22(^5 - ^ 4 ) 

^87 = P \l{L - S5)A 2 i A 2 s s 
^88 = -P\l(L - S $)A2 X A lSf 

Hence 

/ 0 0 0 0 0 0 0 

Z/2| L 22 Z^3 0 0 0 0 0 

L 31 L 32 ^33 1*A 0 0 0 0 

L & ((o) = ^ 41 ^ 42 ^ 43 ^ 44 ^ 4S 0 0 0 

1*1 1*2 1*3 ^54 ^55 0 0 0 

1*1 1*2 1*3 1*A 1-65 1*6 1-61 0 

Lj\ L 72 L 73 In a 1-15 1*6 0 Z /78 

1*2 1* 3 ^84 ^85 ^86 ^87 ^88 

Suppose co is a transmission zero of \|/(ito) : 


then in the notation 


we have 


Suppose for some z : 


\y(io))v = 0 , v * 0 , 

/ 

x ~ b 

A 0 f = to 2 M 0 f 

A b f = B u \ + to 2 M b b ; B%b = 0 

^(to)z - (ti 2 M b L b ((i))z = B„v . 


(L a (co) - to 2 M b L b ((fi))z = 0. 
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Then 


and hence 


which would imply that 


-<D 2 M b b + A b f = 0 
-co 2 A/ 0 / + A 0 f = 0 

co is an undamped structure mode 
B*b * 0 


which is a contradiction. Hence 


where 


/ = jC«D )z 


B*L b (<a)z = 0, 


z * 0 


and conversely. Hence the transmission zeros of \|f(/co) are precisely the zeros of 

Det \B u B*L b ((0)\ = 0 . 

These values of (0 are then the “clamped-clamped” modes of the articulated structure: 

-Oi 2 M 0 f + A 0 f = 0 , /* 0 | 

B*b = 0 . i 

and a subset of these corresponding to b = 0 or L b ((0)z = 0 are the clamped-clamped 
modes of the unloaded structure (every node is clamped). Let us consider first co such that 

Det |L fc (co)| * 0 

so that we can invert L b { co). We have: 

b\ = zi 


b 2 ~ I- 4 \b\ = L22Z2 + L23Z3 


bj — L $\ b \ = L32Z2 + L34Z4 

“ Ltife] = L42Z2 + ^43(23 - Z4) + L45Z5 
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where 


z 5 = ^55(^5 - ^5l z l - ^52 z 2 - ^53 z 3 ~ ^54 z 4) 

These four equations can be solved for z\ , zi, Z3, z 4, in terms of b\ , .... 65, and then 
Z5 can be expressed in terms of b$ , b 4 , £>3 , b-i , b| . Next 

~ Ps 5 ,\ 1 (-^2)^1 l(^S “ s 4)^5 = b 6 = f- i 66 z 6 + ^67 z 7 

67 - ^55,1 1 (^2)^1 l( J 5 “ ^4)^5 = ^7 = ^76 z 6 + f>78 z 8 

- (/’iKL - J 5 )P n (j5 -J 4 ) + / > l 2 (f>-^)/ > 2l(^5 -^4))^5 

= 6 g = Lggze + L^fZi - Lg7Zg . 

These three equations can be solved for z 6 , z 7 , z g : 

LiTL^ibe - ^66 z 6) - ^87^7 8 (^7 - f^76 z 6) *= b% - L%6 2 6 
L%-iL(,-)b(, — L%iL-i%bi — bg = — ^7^78^6 —L%(,)z(, . 

Hence 

Z 6 = (Z^ 7 L 6 7^66 - ^87^7 8^-76 ~^86) '(^87^67 *6 ~ ^87^7 8^7 ~ £g) 

z 7 = ^6l(^6 ~ L 66 z 6 ) 

z 8 = ^7 8 (^7 - ^76 z 6) • 

Hence 

z = L t (w) _1 /!>, Det |L b (0))| * 0 . 

Next let 

**/ = 

Let us determine L* . Now from the form of / 4 fc /, it is convenient to break up L* as 

= ^ 4 , + 
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where 


L A 


2 


-L\b\ 

~L\ ,s 2 ^2 

L\ ,5 j ^3 

L\ jb 4 

0 

~L\ ,Sj ^6 

L 1 , 5,^7 

L\b% 


and {/?,} have been determined in terms of z x . Next let 


L Ai z m h 

so that 

b\ = ~^ 21 z 2 

b-1 = - ^2,S 2 ^2,21 C — 12(' s 2) z l + ^12(^2) z 2] — ^2 ,S S 2 ,2 2 ( _ ^l) z 3 
^3 = ^2,52^52,21 + P\i{ s 2) z t\ + ^2,SjPs lt 2 2 ^1 ) z 3 


^4 = ^ 2 , 5 r {/*$,. , 22 ^t ) z 5 

+ L S T ,2 1 {Lrp)\P li (s T — S 2 )(P 1 1 C^2 ) z 1 + P \l(. s l) z l) 

+ P\l(*T ~ ^l){ P 2\{ s l) z \ + P 22^2) Z 2 + A 2 A 2 ,5 2 ( z 3~ z 4))]} 

h 5 = A 2 ~ + P\iU 2 )z 2 ) 

+ f 3 12(^7 - S 2 ){Pl\{.Sl)Z\ + Pl2{Sl)Z 2 + A 2 'a 2 Si (z 3 - z 4 ))] 

+ ^22(^4 - s t ){P 2 i(s t - s 2 )(P u (s 2 )zi + Pn(s 2 )z 2 ) 

+ Pzz( s T - S 2)(P2\(S2) Z \ + P22{s 2 )z 2 + A 2 ' A 2 _ S 2 ( z 3 ~ z 4)) 

- A 2 1 L l iT [Pii(s t -; s 2 )(Pu(s 2 ) z i + P 12(^2)22) 

+ Pl 2 ^T- S 2 )( P 2 l(^ 2 ) z l + P 22 (.S 2 ) Z 2 + A 2 A 2 S ^Z 3 -Z^))] 
+ >4 2 ^ 2 ,r z 5 }] 
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/l6 = —^2.5 5 {^5 3 ,2 1 (“^2 )( ^11(^5 “ ^ 4)^5 + ^12^5 “ s 4) z 6) + ^5 5 ,2 2 ("*^2) z 7 } 

hi = ^ 2 ,s 3 {^S 5 ,2 1 (^ 2)(^1 \( s S - s A ) t >$ + ^12(^5 “ s 4) z 6 ■+" ^ 5" 5 , 2 2 ( ^2 ) z 8 } 

/lg = ^ 2 {^2l(^ “ Ss)[P\\(Ss ~ S 4)bs + P \l( s 5 “ s a) z 6] 

+ Pll(L - J 5 )[^12(*5 ~ ^4)^5 + ^22^5 “ ^4)^6 + ^2 1 ^2,5 5 ( z 7 ~ Zg)]} • 

Hence finally, in terms of L A and L bi (2.2) becomes 

{-<S?M b L b + L a + Yi(0B M B*4)z = B tt v . 

For 0) such that L b (<o) is nonsingular we can write: 

u - B*b = B*L b (-<Si 2 M b L b + L A + yi <*B u B* u L b y X B u v 
= B*(-wX + L A L b l + yi<oB u B*)- X B u v . 

Here 

B u is 48 x 8 
£* is 8 x 48 

and 

(~G) 2 M fc + L A L b X + Y/toS u fi*) _1 is 48 x 48 
and is conveniently broken up into 6x6 blocks, denoted 

D-{Dij}, i.j- 1 . .8 

We can now calculate 
Now 

B u 

B u u = 

B u 

where each B u . is 6x8, and letting 
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“1 

“2 

«8 


we have 


B u ^u = col [0, u u u ly 0, 0, 0 ] 

B Ul = = 0 

B Ut u = col [ u 3 , u 4 , 0, 0, 0, 0 ] 

B Ui u = col [ 0, u 5 , « 6 . 0, 0, 0 ] 

- 0 

B u u = col [ « 7 , u 8 , 0, -40«s, 40 UT, 0 ] . 

8 


Hence letting 




b = 


*>8 





+ B* b A + B* 




8 


and writing 


u 

V 

w 



$2 

<t>3 


Z^/i = col [ v, w, 0, 0, 0, 0, 0, 0 ] 

B* h = col [ 0, 0, u, v, 0, 0, 0, 0 ] 

w 4 

B* } h = col [ 0, 0, 0, 0, v, w, 0, 0 ] 

B* h = col [ 0, 0, 0, 0, 0, 0, u + 40<{>2 > v-40<5>, ] . 

u 8 
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Hence 


(B*DB U )V = B* f E D t jB Uj v + lD 4y fi„.v + B* } YD 5J B Uj 


= B u x + ^14 B u f + B> 15 B u j + £> I8# U ,) V 

+ B* t (D 4l B Ui + D 44 B Ut + D 4 $B Ui + D 4 gB Uf ) 

+ B u i ( D n B u i + d 54® U4 + 0 55 B Bj + D 5i B Ut ) 

+ B Z,{Dsi B Ul + D&4 B Ui + D iS B u s + °88®u,) 

This shows in particular that we do not need to calculate all the 
at location S it i = 1, 4, 5, 8 are given by 


B * Ui (Dn B u, + D iA B Ui + D iS B U} + D ig B B| )v = fi u *.D M B u .v + 


where the first term involves only the sensors at locations 5, 
represents the coupling to sensors at other locations. Also 


B ZD n B u: 


D, | 22 v 1 + I ,2 3 v 2 

D 1 1 ,3 2 v l + D \ 1 ,33 v 2 
0 

0 


V = 


where 

^ii = {Dn.ij } > 


Similarly 


B t t D^B Ui 


0 

0 

^44,1 1 v 3 + £> 44,I2 v 4 
£>44,21 v 3 + £) 44,22 v 4 
0 
0 
0 
0 


ij = 


8 

v + B* t ’LD Sj B Uj V 


. Also the controls 



and the summation 

vi 

v 2 

Vg 
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0 
0 
0 
0 

^55,22 v 5 + ^55,23 v 6 
^55,32 v 5 + ^55,33 v 6 
0 
0 

0 
0 
0 
0 
0 
0 

(^88,1 1 + 40O 8 8 P 1 5 ) v 7 + (^ 88,1 2 - 40D gg , 4 )v g 

+ 40{(D gg _ 5] +40D g8i55 )v7 + (£ > gg 52 -40D gg 54 )v g } 

(£*88,21 +40D g 8 25 ) v 7 + (£*88 ,22 8 ,24 ) v 8 

+ 40{(D gg 41 +40D gg(45 )v7 + (£ > 88>42 -40Dg g 44 )v g } 

It must be noted that in terms of feedback 

v, = v( 0 ) 

v 2 = W( 0 ) 

V 3 = U(S T ) + 100$ 2 (5t-) 
v 4 = v(S T ) - 100<j>j (S T ) 
v 5 = v(S 4 ) 

v 6 = w(S 4 ) 

V 7 = it(L) + 40$ 2 (L) 

v 8 = \(L) - 40<j), ( L ) , 

the dot denoting derviative. 
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Let 


7(/co) = L A (a>)L b (a>) 

so that the poles of T(i(o) are the zeros of L b ( (0). We shall show that T(i(0) is self-adjoint 
but not nonnegative definite! Let b\ , bi be arbitrary real, and let 

L fr (co)'' 6, = Zi 

4(G))" 1 4 = z 2 

/, = £(o ))z, 


Then 


h = £(«) z 2 


X\ = 


/l 


*2 = 


/ 2 
4 


= [Aq/,,/ 2] + [A b /, , 4] = G) 2 [Mo/,,/ 2 1 + 14/, ,41. 


Similarly 


Since 


[Ax 2 ,x l ] = [4/ 2 ,41 + 0 ) 2 [M 0 / 2 ,/,]. 


[AjCi,X 2 ] = [Ax 2 , JC! ] 

we have 

l4/ 2 . 41 = (4/l , b 2 ] 

or, 

[4/i.4J = [^(a)) 4 (( 0 )" , 6 1 , 4 ] = [ 4 / 2 . *> 1 ] = [4(®)4(®)'^2. 4]. 

Hence T(i(x>) is self-adjoint. Next 

[T(i(£>)b x , by] = t4/l» 41 = [ Ax x ,x x ] - G) 2 l M 0 /, , . 

Hence 7(/C0) is not nonnegative definite. Note finally that the structure modes frequencies 
can be expressed: 

Det |-G i 2 M b + T(/G))| = 0 . 
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